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Abstract 



The contribution of subleading jet functions to inclusive decay distributions of B 
mesons are derived from a systematic two-step matching of QCD current correlators 
onto soft collinear and heavy quark effective theory. Focusing on the tree level matching 
of QCD onto soft collinear effective theory, the subleading jet functions are defined to all 
orders in as{ni) (with ~ mj,AQCD) and are calculated explicitly at first order in as{ni)- 
We present explicit expressions for the decay rates of 5 ^ I v and the Qy-.^ — Qi^^ 
contribution to B ^ Xg 7, where the subleading jet functions are multiplied by a tree 
level hard function and appear in a convolution with the leading order shape function. 
Together with the recent two loop calculation of the leading order hard function for 
B — > Xu I this paper will allow for a more precise description of inclusive B decays in 
the end point region. 



1 Introduction 



Charmless inclusive B decays play an important role in our understanding of the standard 
model and its possible extensions. The inclusive semileptonic B decay currently allows for 
the most accurate determination of \Vub\, one of the fundamental parameters of the standard 
model, while the inclusive B ^ Xg^ rate is used extensively in constraining models of new 
physics. 

Since Aqcd is much smaller than the 6-quark mass {rrib), one would expect that various 
physical observables for inclusive B decays can be expressed in terms of a local operator 
product expansion (OPE), where the various operators are suppressed by an increasing power 
of nib. This is the case for the partial and total rate of B X^lv, where schematically we 
have 

i=2 j b 

The current state of the art is that Cq is known at 0{a'^) [1, 2], while Cg [3] and [4] are 
known at C(q;°). For C2, the coefficient of the "kinetic energy" operator c^, is known at 0{as) 
[5], while c|, the coefficient of the "chromomagnetic" operator, is known only at 0{a^g) [6, 7]. 

For the charmless inclusive B decays, S — > X^. 7 and B — >■ Xu 1 1', the situation is more 
complicated. Experimental cuts force the hadronic jet X to have large energy Ex ~ nib, but 
only moderate invariant mass Pj- ~ ni^ Aqqb- Consequently there are three energy scales in 
the problem: a hard scale {iih ~ ^b) , a hard-coUinear scale ~ ^/mbA^qcD) and a soft scale 
{lis ~ Aqcd)- For this kinematical region, often called the "end point region" or "the shape 
function region", the partial rates can be expressed in terms of a non-local OPE. Thus, for 
B Xu I V and the Qi^ — Qi^ contribution to S ^ X^ 7 we have an expansion which is 
schematically [8]: 



i=l ^b 




where the hard functions (if, hi) and the jet functions (J, jj) are calculable in perturbation 
theory, while the shape functions (5, Sj) are non-local light-cone operators which are non 
perturbative objects^. Factorization theorems such as (2) are most conveniently proven using 
the Soft CoUinear Effective Theory (SCET) [12, 13, 14]. 

The current state of the art is as follows. if„ was recently calculated at 0{a^^ [15, 16, 17, 18] 
and Hg [19] is known^ at 0{a^. The leading order jet function J is known at 0{o?^ [25]. Of 
the 1/mb corrections, only the terms of the form Yb{ ■ j\ s\ are explicitly known at 0{c^^ 

^For other contributions to rfFs, such as Q^-^ — Qf,g, one finds that more complicated factorization theorems 
hold [9, 10, 11]. We will not discuss these contributions in this paper. 

^The 0(q:^) expression can be extracted from known results in the literature [20, 21, 22, 23], as was done 
in [24] for a "normalized" at the scale [ih = mb- 
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[26, 27, 28]. Therefore the factorization formula was proven only for the leading order term 
and one of the l/m^ suppressed terms. 

The knowledge of the leading order hard and jet functions at 0{as) and the h'l ■ (E) s\ 
terms at 0{a'^) (as well as known, but not properly factorized, as/ ml and l/ml corrections) 
was the basis of the precision determination of \Vub\ in [29, 30, 31]. In order to improve 
the accuracy even further, one would like to know as much as possible about the properly 
factorized a^, as/rrih, and the corrections, in decreasing order of importance. What 

would the calculation of these corrections entail? 

In order to find the least important term, namely the l/ml corrections, the heavy-to- 
hght SCET currents need to be matched at tree level to fourth order in the SCET expansion 
parameter KQCTt/mi,. The most important term, namely the corrections, should appear 
shortly now that the last ingredient, Ha, was calculated at 0{a1) [15, 16, 17, 18]. The second 
most important correction is of order ag/mi,. From (2) one would naively expect to find three 
terms which scale like l/m^ and need to be calculated to 0{as): h} ■ ® s", ■ ® s", and 
/i° • j° (8) (the subscript 1 is implicit). Less formally, we expect to find subleading hard, jet 
and shape functions, respectively. Let us discuss each of these terms. 

Since the hard functions are products of Wilson coefficients extracted in the matching of 
QCD onto SCET, they depend only on kinematical quantities which scale like m^, i.e. mj, and 
Ex + \Px\- As such they always scale as C(l) in the l/mb expansion, so it is clear that terms 
of the form ■ j° (g) s° cannot appear at any order in a^. In other words, subleading hard 
functions can appear only when they are multiplied by subleading jet or shape functions. 

The "subleading shape functions" (SSF), i.e. terms of the form /i" ■ j° (g) which arise 
already at 0{a'^), can be calculated, in principle, at 0{as), i.e. both and j*^ need to be 
calculated at 0{as). For the former a one loop matching of the SCET current to second order 
is needed, while the calculation of the latter was outlined (but not exphcitly done!) in [26] 
and [28]. 

The focus of this paper is to prove that terms of the form ■ ® <s°, i.e. "subleading jet 
functions" (SJF), indeed exist in the factorization formula for inclusive B decays. We will first 
prove the existence of such terms by showing that the partonic 0{as) terms in the hadronic 
tensor which are l/mb suppressed in the end point region arise from two momentum regions: 
soft and hard-coUinear. We will then show how the soft region is accounted for by the parton 
level one loop diagrams of the known O(a^) h'^ ■ j° ® term, and reproduce the hard-collinear 
region via time order products (TOPs) of subleading SCET currents. After establishing the 
need for the /i" • (8) term, we will calculate the subleading jet functions via the usual two 
step matching. In the first step the QCD currents and Lagrangian are matched onto SCET at 
tree level and to second order in \/A^cd7^- the second step the SCET current correlator 
is matched onto Heavy Quark Effective Theory (HQET) [32] and the subleading jet functions 
are extracted. Since we are interested in as/mb suppressed terms, and the subleading jet 
functions start at C(q;s), it is sufficient to consider only the case of tree level matching of 
QCD onto SCET, for which we can use known results from the hterature. For this case we 
will define the subleading jet functions to all orders in ^^(//j), and calculate them explicitly at 
first order in Q;s(/ij). 

The subleading jet functions' contribution is, in some sense, the most important term at 
order as/mb- When integrating over larger portions of phase space, the one loop subleading jet 
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functions' contribution is no longer power suppressed. The other as/mi, term in the factoriza- 
tion formula, namely the 0{as) ■ (g) contribution, although formally as/rrii, suppressed 
in the end point region, is expected to become even more power suppressed and thus is less 
important outside of the end point region. In other words, the terms which we will calculate 
are kinematically and not hadronically suppressed, and as such are important outside of the 
end point region. Incidentally, experiments are starting to probe the kinematic area outside of 
the end point region. Finally, since the subleading jet functions appear in convolution with the 
leading order shape function, their inclusion does not introduce new hadronic uncertainties. 

Another motivation for studying subleading jet functions is that they arise also outside of 
the context of flavor physics. Since these functions encode the interaction of hard-coUinear 
quarks and gluons, the same functions are expected to appear in the x — > 1 region of deep 
inelastic scattering [33, 34, 35, 36]. 

The rest of the paper is organized as follows. After a short review of known results in 
section 2, we calculate in section 3 the partonic 0{as) terms in the hadronic tensor which are 
l/rrib suppressed in the end point region. We then show that they arise from two momentum 
regions: soft and hard-coUinear. In section 4 we explain how the soft region is accounted 
for by the parton level one loop diagrams of the known ■ j° terms, and reproduce 
the hard-collinear region via time order products of subleading SCET currents. In the main 
section of the paper, section 5, we define the subleading jet functions, for the case of a tree 
level hard function, to all orders in as{fJ.i) and calculate their one loop expressions. After a 
short discussion of their renormalization, we present properly factorized expressions for the 
decay rates of S — > X^lV and the Q'j-y — Qj-y contribution to B -^s7- In section 6 we 
present our conclusions. The appendices contain proofs for some of the statements made in 
section 5. A reader who is mostly interested in the phenomenological results, can skip section 
3 and 4 and proceed directly to section 5. 

2 Review 

In order to make the paper self-contained, we review in this section some known results 
about inclusive B decays, as well as some basic ingredients of SCET. For a more detailed 
account see [37]. 

Kinematical Variables The kinematics of inclusive B decays is such that in its rest frame, 
the B meson decays into a hadronic jet carrying momentum Px and a non-hadronic jet (a 
lepton pair for B — > 1 9 and a photon for B ^ Xg^) carrying momentum q. Denoting by 
Mb the mass of the B meson and by v its four- velocity, we have therefore M = Px + q. 
Taking the four velocity of the B meson to be v = (1, 0, 0, 0) and q to point in the negative 
z direction, we define two light-like vectors = (1,0,0,1), = (1,0,0,-1), such that 
n + n — 2v, n ■ n — 2, and n • v — n • v — 1. Any four vector can be decomposed as 

— n ■ a— + n ■ a— + a']^. (3) 
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Notice that we have taken v± — q± — 0. Having fixed the two hght-hke vectors, rotational 
invariance imphes that the transverse indices can only be contracted using 

9T = r- \ , = \e^^-'n^n,, (4) 

where 69123 = 1- 

Conservation of 4-momentum implies that for the i? ^ 7 decay mode there is one 
independent kinematical variable, which we can take to be the photon energy or n ■ Px ~ 
Mb — 2E^. The B ^ X^lV decay mode has three independent variables which we can take 
to be [29, 38, 39] 

P+ = Ex-\Px\=n-Px, P-=Ex + \Px\=n-Px. Pi = Mb - 2Ei. (5) 

These are the "hadronic" variables. It is also useful to define a "partonic" set of variables. 
Let A = Mb — rri},, where is the b quark mass. Defining p = rribV — q = Px — Av, we have 
n-p — n-P — A, n-p — n-P — Ae^ the corresponding partonic variables to P+ and P_. 
The hadronic tensor is naturally expressed in terms of n • p and n ■ p. Notice also that by 
construction p± — 0, and as a result p^ — n ■ p n ■ p. 

Hadronic Tensor Partial rates for the inclusive decays B — > I i^, and the Q^^ — Qr^ 
contribution to S — > 7 can be calculated using the optical theorem. The central object to 
consider is the hadronic tensor, which is the discontinuity of a forward matrix element of a 
correlator of two currents: 



n2MB 



lm{B{v)\i J rf^xe^«-^T|j/(0) Jj(x)} (6) 

where, again, v is the velocity of the B-mcson and q is the momentum of the lepton pair 
(photon) in the B X„ lu {B ^ 7) decay. The currents can generally be written as 
j/ — bViq and Jj = qTj b. For semileptonic decays Fj = — 75) and Tj = 7*^(1 — 75), and 
for the — Qr-y contribution to S ^ X^ 7, Fj = |(1 — 75)7^ fl, F^ = |fl7^(l — 75). In order 
to somewhat simplify the traces in the expression for the hadronic tensor, we assume that Fj 
and Tj contain the same number of Dirac's gamma matrices, but otherwise we take Fj j to be 
arbitrary Dirac structures. 

For radiative decays the hadronic tensor is given in term of one function W = Wij. The 
Qi-y — Q77 contribution to photon spectrum can then be written as^ 

^i?£^?IWrm^|C??riy(P^). (7) 



dE^ Att^ 

For semileptonic decays we can decompose the hadronic tensor in terms of five functions, 

-(rW2 + v^'v^Wa + {n^v" + rfv^)Wi + n^n^'W^. (8) 



^See [29] for the exact definition of the various parameters in the this equation. Notice that W equals 
-2U{iih,tii)r^ of [29]. 
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The triple differential decay rate can be written in terms of W^i, . . . , W5 as [29] 



_ G'jp\Vub\ 
dP+ dP^ dPi ~ 167r3 



{Mb - P+) 



(P_ - Pi){Mb -P- + Pi- P+) W^ 



+ {Mb - P-){P- - P+) ^ + {P- - Pi){Pi - P+) + W,+ -W, 

2 \4: y 



(9) 



where 



P^-P. 



y = 



+ 



Mn-P^ 



(10) 



Known l/m^ corrections The hadronic tensor can be factorized as in equation (2). In 
this paper we will be interested in the terms which are suppressed by one power of m^. There 
are currently two types of these terms which are known. The first type are "subleading shape 
functions" i.e. properly factorized terms of the form • s^, where and j'^ are explicitly 
known at The second type are "kinematical power corrections", i.e. terms calculated 

within the parton model which arc suppressed both by and l/mj. These terms will be 
properly factorized in this paper. We now briefly review these two types. 

The contributions of the form /i° • j° (8) to the hadronic tensor, i.e. the subleading 
shape functions, were calculated using SCET in [26, 27, 28] (for earlier partial calculations see 
[40, 41, 42, 43, 44]). Here we use the results of [27]. The above contribution to the hadronic 
tensor is: 



W 



- SSF 



du 6{n ■ p + uj) 



duj 6{n ■ p + u) 



mi, nib 

fu{^) rj. , fv{uj) ' 

— -t 1 + — -1-4 

n ■ p n ■ p 



n ■ p 



n • p 



viu) 
n ■ p 



where 



tr 



To = - tr 



1 + f 



tr 



Ta = - tr 



ri^75r,-i^(y-yi 
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(12) 
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The subleading shape functions are defined as: 



and 



mb{i J (iST{/i(0) [0,a;_]/i(a;_)£f (^)}) = J rfa;e-^'"^-^s 



n-x/2 





n-x/2 



/ dw e~ 


■^'""■^ S{u) , 


j dw e~ 




f 




j duj e~ 




I dw e~ 





n-x/2 n-x/2 



ti 

= J due-^-'^-^Uiuj), 

n-x/2 n-x/2 



ti 

= / duje-i^'^-^f.iuj), (14) 



where k, I are color indices, 5" in equation (14) is a soft Wilson line defined in [27] (not 
to be confused with the leading order shape function S{iju)\), [x,y] = S{x)S^{y), JC^^^ is the 
next-to-leading term in the expansion of the HQET Lagrangian, and 

^ ^ {B{v)\h...h\B{v)) 

The second type of terms, namely the kincmatical power corrections, can be found in [29], 
where they were called In that paper the corrections were convoluted with the "tree 

level shape function" , in absence of proper factorization. The relevant expressions are, for the 
Qv-y ~ Qv-y contribution to S — > 7, 

W = -M^p^^^^ dwS{u,y^,) (-15 - 16 Inx) , (15) 
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and for B — > I ly, 



^km(l) 



^kin(l) 



71^3 + W4 + -W5 

4 y 



kin(l) 



1 






Mb- 


P+ 


47r 


1 




Cpasift) 


Mb- 






1 







where 



du S{u}, ^i) 



duj S{u}, Hi) 



y \y ) ^. 



\ 22 8, y 
4 + - In - 

y y X. 



(16) 



X — 



The "hatted" function S{u,iJ,i) is related to S{lj), defined in equation (13), by a change of 
variables and a 1/mf, suppressed term. The exact relation can be found in [29]. 

Since the expansion in [29] was organized in inverse powers of Mb — P+ instead of nib, for 
future reference we will need also the leading order term for Wi 



where 



Hui{y,fJ'h) = H- 



rP+ 

{y,l^h) / ( 
Jo 



-4 In 



du; yrub J{ymb{P+ - cu), Hi) S{u;, Hi) , 



2y^b , -t^.^^yrrib 21n|/ _ , tt 



(17) 



+ 101n- 



41ny 



AL2il-y) 12 

l-y ^' & 



and 



J(p^/.)=5(p2) 



+ 



At: 



4 In 1^-3 



(18) 
(19) 



Some Ingredients of SCET SCET is the appropriate effective field theory to discuss inclu- 
sive B decays in the end point region. The SCET expansion parameter is y/X = a/Aqcd/^- 
For the following we will need the expansion of the SCET heavy-to-light currents to second 
order, as well as the leading order hard-collinear Lagrangian*^. These are most conveniently 
hsted in [27]. 

First, the hard-coUinear Lagrangian is 



(0) 



in ■ D 



_ (0) 

(0) 'r±hc 

he 



■ip 



(0) 



(20) 



^The power corrections to the hard-coUincar Lagrangian always involve extra soft particles (apart from the 
heavy quark) and as such do not contribute to the subleading jet functions. 
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where the Lagrangian is written in terms of "sterile" fields, i.e. after the soft degrees of freedom 



were decoupled via a field redefinition. In the last equation iD^^^'^ — id^ + gA^P}^ and 



W^Pexj^iig j dtn ■ Afl{x + tn)\ 

\ — rv-i ' 



(21) 



Next we need the expressions for the currents. We present them in terms of the "calligraphic 
fields" [45] 

X = lyt^(o) ^ y^M^ ^ [iDfl^W) . (22) 

Generally speaking, in the SCET expansion of the currents (and the Lagrangian), the 
power suppression arises from two separate sources, not mutually exclusive: presence of power 
suppressed components of the hard-collinear gluon field or the hard-collinear covariant deriva- 
tive, and presence of soft fields or their covariant derivatives. For the purpose of this paper 
we need only currents of the first type, which are 



'I 

i 



2 in - d 2m6 



J(2) 



— X V— — n ■ Ahc [S'^h) — X r — — - n ■ Ahc (S^h) 
zrub in ■ a ^ 



xr 



in ■ d 



{if±hc 4±hc)_ ^^t/j^) _^ ^ '^^-L'^c 



^-^l;^7;4Ms'h)^ . (23) 

in ■ d ^ 



where we have suppressed the overall e ^^bv-x f^c^or in each term. For completeness we list 
also the currents of the second type, which we will not use 



J 



t{2) -y p 

■^not used — ^ 



n ■ X 



X , X 



iS^n.Dh)^_ + --^iS^D,D^h)^_+ 



2mh 



h 



+ x^i^^^_rx'i{s''D,h)^ -x(^4^uc^^r + Tj-4^^)jx^^{s^D,h)^ . 

^ in ■ d \ ^ in ■ d / 



The second term in J, 



(1) 



not used 



(24) 



does not contain any soft fields apart from the heavy quark. Still, 



we can ignore its contribution in this paper, since we can always set p±, where p is the total 
hard-collinear momentum, to zero. 

It should be noted that these currents were matched at zeroth order in as{^h) and as a 
result the Wilson coefficients are always equal to 1. When matching beyond zeroth order 
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in as{iJ,h), one would expect more complicated currents with multiple non-localities, see for 

example the one loop matching onto the first order SCET current in [46]. The resulting 
contributions to the hadronic tensor would be suppressed by as{^h) x Oisil^i) X ^/fT^b and as 
such are much smaller than the contributions considered in this paper. 



3 Analysis by regions 

In this section we will calculate, within the parton model, the one loop corrections to the 
hadronic tensor, defined in equation (6), that scale as 0{X^) in the shape-function region, 
i.e. terms which are suppressed by l/m^ compared to the leading order terms, which scale as 
0{X~^). We perform the calculation for a general Dirac structure in the hadronic tensor. In 
order to simplify the calculation we neglect the residual momentum of the b quark, i.e. we 
work with on-shell b quarks. As a result, all the terms are constants or logarithms of the form 
In (n ■ p)/{n ■ p) = In r, where p is the partonic momentum of the jet, defined in section 2, 
which satisfies p± = and p^ = n ■ p n ■ p. 

We perform the calculation both in "full QCD" and by using the method or regions [47, 48] . 
We find that only two kinematical regions are needed: a hard-coUinear region, where the loop 
momentum scales as (1, A, A^/^), and a soft region, where the loop momentum scales as (A, A, A). 
There is no contribution from a hard region, were the loop momentum scales as (1,1,1), which 
is in accordance with the lack of terms of the form ■ j° s° in the factorization formula. As 
explained in the introduction this is a general result which holds also beyond one loop order. 
Typically, we find that terms of the form In r can be decomposed as: 



In r = In 



= — -Mn^ + -+ln ^ — ^ 
n ■ p e p'^ e jj,"^ 



hard— coUinear 



(25) 



At one loop there are several diagrams that contribute to the hadronic tensor. These diagrams 
are shown in figure 1. We present the results for each diagram separately, namely, the "Self 
energy" diagram (top left), the "Box" diagram (top right), and the "Vertex" contribution 
which is the sum of the diagrams on the bottom line of figure 1. We use Feynman gauge 
throughout this paper. 



3.1 Full QCD 

Calculating the diagrams we find the following results: 
"Self energy" : 



W,,^^e{p')tr 



1 1 



n • p 4' 



(26) 
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Figure 1: One loop diagrams contributing to the hadronic tensor, top left: "Self energy" 
diagram, top right: "Box" diagram, bottom line: two "Vertex" diagrams. The letter next to 
each solid hne denotes the flavor of the quark. 



"Box" 



w^, = '^e{p'){ tr 



47r 



n • p 



;-l -Inr) +tr [Fi^Tj 



1 1 
nib 4 



11 Tl ' U \ 

+ tr [r, r, ^] — - ( -2 - In r) + tr [F, ?^ r, ?i] - , (27) 



"Vertex" 



tr 



^%2)ar[F,^F,.^] ^l(r + ln 



1 1 /3 



F,^7^7^F,-7^7^ 



rUb 4 V2 

1 + , 

+ tr 



tr 



Fi «i F, 



1 + / 



n ■ p 



;i + 21nr) 



r.727^?lr,i±^7/3 7^ 



1 1 

16 



tr [F,7^F,7^] 



1 1 
rub 8 



(28) 



In order to check these results, we can compare them to the expansion of the one loop expres- 
sions of the hadronic tensor for B — > I i^, and the Qr^ — Q-j^ contribution io B ^ Xg 7. 
This is most easily done by using equations (15)-(19) in section 2, which are taken from [29]. 
In that paper the correction were convoluted with the "tree level shape function". We can 
undo this convolution by the replacements: 



\ rrib, y 



n ■ p 
nib ' 



X 

y 



n ■ p 

n-p' Jo 



duj S{u, Hi) — > 1. 



(29) 
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For B — > Xulu we also need to expand Wi*^ in powers oi n ■ p/n ■ p. In total we find for 

^^"^ 2 (-15 - 16 In r), 



and for B ^ Xulv. 



W 



Wi 



Wo = 



47r mft 



Cpas / 10 9 12 ^ 4 ^ 

— in r H in r 

An \mb n-p n-p nib 



(30) 



4 y 



Cpas 2 
Att n-p 

Cpa, / 4 _ 22 



Inr 



. . _ _ (31) 

47r ymj, n-p n-p / 

Summing over (26), (27), and (28), and calculating the traces for each decay mode using the 
expressions after equation (6), we find complete agreement with (30) and (31). 

We are now ready to repeat this calculation using the method of regions. We perform the 
calculation in d = 4 — 2e dimensions and use dimensional regularization to regularize both the 
UV and IR divergences. We also implicitly take /j, — > iJ,e""^^'^{4:7i)~^^'^. 



3.2 Hard-collinear region 

For the hard-collinear region we find the following results. 
"Self energy" : 



W,^eip^)^tr 



i + i 



1 1 



n ■ p A' 



"Box" 



Cf^s 

At: 



tr 



n ■ p \ e 



(32) 



+ tr [P. ^ r^. ^] i- 1 (^-1 - 1 - In ^) + tr [F, ^ r, ^] ^ 1 ^ , (33) 



"Vertex" 



tr 



F,:^F 



1 1 n\ 
n ■ p \ e 2 p"^ 



+ tr [F, ^ F, fl] i- i Q + ^ + In ^) - tr [F. 7! T, 



1 1 

nib 8 



- tr 



F,^7^7^F,-7^7^. 



1 + ' 



+ tr 



1 1 

nib 16 



(34) 
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3.3 Soft Region 

For the soft region we find that the "self energy" diagram does not contribute. For the other 
diagrams we have: 

"Box" 



Wi, = 9{n ■ p) 



tr 



n ■ p \e 



-3-ln 



+ tr[riyirjfi] — ^ + tr[ri^rj-^] — ^ ( - - 1 -in 




(35) 



"Vertex" 



Wij^9{n-p)^^<tT 



1 + 



n ■ p 



2 - - 1 - In 



(n ■ pY" 



+ 



+ tr r,^r,-^ — - — + ln^^ 



(36) 



Adding up the two types of regions wc find that, as expected, the sum of the hard-collincar 
and soft regions reproduce the full QCD result. Notice also that the structure of the soft 
region is simpler than that of the hard-coUinear region. In the next section we will see that 
the reason is that the soft region is accounted for by the parton level one loop expressions 
for only two subleading shape function, while for the hard-coUinear region we need several 
subleading jet functions. 



4 Effective field theory calculation 

We have seen in the previous section that the one loop corrections to the hadronic tensor which 
scale as 0(A°) in the end point region, arise from two kinematical regions: a hard-coUinear 
region and a soft region. In this section we will see how the soft region is accounted for by 
the one loop corrections to the subleading shape functions' contribution calculated within 
the parton model. The hard-coUinear region is accounted for by the contribution of the time 
ordered product of power suppressed SCET currents. This calculation is the basis for the 
subleading jet function calculation which we perform in the next section. 



4.1 Soft contribution 

The contribution of the soft region can be fully accounted for by calculating, within the 
parton model, the one loop corrections to the subleading shape functions. In particular there 

is no need to introduce new subleading shape functions. From equation (11) we see that the 
hadronic tensor depends on several subleading shape functions. Naively, one would assume 
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Figure 2: One loop diagrams contributing to the parton level expression of u{u!). 



that we need to calculate the one loop corrections for ujS{uj), s{uj),t{uj),u{uj),v{uj), as well as 
the four-quark shape functions /^(cu). In practice, only the contributions of ci;5'(a') and 

u{u!) are needed for the following reasons. 

- We have chosen the coordinate system such that v± = 0. As a result the matrix elements 
of the operator corresponding to t{uj) and v{u), vanish at one loop, since they only 
contain gluons which have perpendicular polarization. 

- The matrix elements of the operators corresponding to /„(c^) and /^(a;) vanish since 
they involve scaleless integrals over the n component of the light quark momentum. 

- Setting the residual momentum of the heavy quarks to zero, we find that the matrix 
element of the operator corresponding to s{lj) vanishes. 

As a result we need the "zero external gluon" matrix elements of the operators corresponding 
to ljS{uj) and u{lj). The first can be extracted from [38]. After setting the residual momentum 
k to zero, we find 



UJ 



TT 



— + ln— + 1 



(37) 



For u{u!), we need to calculate the one loop amplitude which is the sum of the diagrams 
shown in figure 2. The relevant Feynman rules needed for this calculation involve zero and 
one external gluon. We have 




—k']_ d{(jj — n ■ k) 



+ 



n ■ I 

fa g^fj. 

(n • ly 



{k - 1)1 5'{ou -n-k + n-l)-kl 5'{ou - n ■ k)j 



{2k± • /± - /i) S{uj-n-k + n-l)-S{uj-n-k) . (38) 



The notation is such that k is the incoming heavy quark residual momentum and /, /x, a are the 
outgoing gluon's momentum, polarization, and color index, respectively. For the one gluon 
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Feynman rule we have omitted terms in which the gluon has a perpendicular polarization, 
since such terms do not contribute to the one loop amplitude for Vx_ = 0. 

These Feynman rules were first calculated by Trott and Williamson in [49], where the 
corresponding operator is called QsIlJjT). Keeping only the terms which are proportional to 

in the one external gluon Feynman rule of [49], our result agrees with theirs accounting 
for an overall factor of 2/m{, and a different sign of the heavy quark residual momentum, 
both arising from the slightly different definition of u{uj). Calculating the one loop amplitude, 
which was not calculated in [49], and setting the residual momentum k to zero, we find 



u 



parton 
bare 



TT 



(39) 



Inserting (37) and (39) into (11), and setting s(a;), t(a;), v{uj), fuij^) and /^(o;) to zero, we 
find: 



W^f^e{n.p)^ < tr 



1 ^3_3^(^_5 



n ■ p \e 



.,[r,^r,w-ffli-i(-Uiu<!i^ + i 



(40) 



which is the total contribution of the soft region, i.e. the sum of equations (35) and (36). 

At this point we should note that the question of operator mixing and renormalization with 
regard to u{uj) is still open, as it was not considered in [49]. The main complication being the 
need to introduce new subleading shape functions which u{uj) can mix into, and establishing 
the closure of the basis. The analysis of this question goes beyond the scope of this paper. For 
our purposes the important point is that the terms of the form ■ j° ® s^, that were already 
calculated in the literature, reproduce the contribution of the soft region. As a result in the 
final factorization formula, no new terms are needed to account for this contribution. 



4.2 Hard-CoUinear contribution 

In order to reproduce the result of the hard-coUinear region, we need to consider various 
combinations of the SCET currents and the SCET Lagrangian^. Symbolically we need the 
following combinations: 

Jt(l) J(l)^ Jt(2) J(0) + jt(0) j(2)^ jm j(0) J ^4^^(0) ^ _^t(0) J(2) J ^4^^(0)^ ^^^^ 

where C^^^ and J*^*) are defined in equations (20) and (23), respectively. Calculating the 
contribution of each combination we find that the non zero combinations are: 

^It is easy to show that the TOP of a first order current with the zeroth order current vanish since the 
perpendicular components of the hard-coUinear momentum can be chosen to be zero. This corresponds of 
course to the absence of 1/ ^/fnb corrections to the hadronic tensor. 
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Jt(l) J(l) 



+ tr 



An 



tr 



i + f 



^ 2 

1 + V' 



+ tr [Fj ^ ^] ^ ^ 



n • p 4 

+ tr 



ml 16 



1 1 
rrib 32 



jt(2) j(o) + jt(o) j(2)^ j(2) ^ _xrJ^n-Ahc{S^h) 



1-ln^ 



jm j(o) + jt(o) j(2)^ ji2) ^ _xrj^n-Ahc {S^h}^ 



Cpas 
An 



tr 



Fi^F, 



1 + f 



^ — + 2-ln^ 



Jt(2) J(0) J ^4^ ^(0) ^ jt(0) J(2) J ^4^ ^(0) ^ j(2) ^ _ ^ V ^ U ■ Aha {S^ h) 



tr [Fj yj, Fj yi] 



11/1 3 



In 



Jt(2) J(0) J ^4^ ^(0) ^ jt(0) J(2) J ^4^ ^(0) ^ j(2) ^ _ J F ^ 71 • ^1^, (^t/^) ^ 



4n 



tr 



n ■ p \ e 2 



— I In — 



p^ 



jm j(o) J ^(0) ^ jt(o) j(2) J ^(0) ^ j(2) ^ _ X r ^ {ii±Hc4^Hc) ^^t/,) 



An 



tr 



- tr 



F,^7f7^F,7^7^ 



1+i 
' 2 

1+y 



\ e 2 



+ tr 



1 + 



1 1 

nib 16 
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Jt(2) J(0) / d'^Z Cf + Jt(0) J(2) / ^(0) ^ j(2) = X %^ ^ I r I [S^ 



An 



tr 



1 + i 



+ tr 




Using 7^ = 7^ — n^ifi,/2 — n^i/i,/2 and 2ili = iji + ifc^ we can show that the sum of (42) — (48) is 
equal to the hard-colhnear contribution i.e. the sum of (32), (33), and (34). Again we see that 
the contribution of the hard-colhnear region is more complicated than that of the soft region. 
As we will see in the next section there are seven different jet functions that contribute at one 
loop, compared to only two subleading shape function that are needed to reproduce the soft 
region. 



5 Subleading jet functions 

Following the calculations of the previous sections, it is clear that in order to properly factorize 
all the terms in the hadronic tensor that are both as and l/rrib suppressed in the end point 
region, we have keep the subleading shape functions' contribution and replace the so called 
"l/mj kinematical corrections" of equations (15) and (16) by the contribution of the subleading 
jet functions. By "jet function" we mean the discontinuity of a Fourier transform of a vacuum 
expectation value of a time ordered product of hard-collinear fields. By "subleading" we mean 
that these functions scale as C(A°) in the end point region. 

Before going into the details of the analysis of each subleading jet function, we wish to 
make some general remarks. The first issue is factorization. In the following we establish 
factorization formula for terms of the form ■ (g) s^. More explicitly, we always have = 1 
and s° is the leading order shape function. The hadronic tensor factorizes as 

Wij = J] C„ tr [F, ... F,- ...] f dcujaipl) S{uj). (49) 

Here Ca is a simple kinematical factor of mass dimension -1, e.g. 1/n ■ p, 1/mb and the 
argument of the jet function Ja is = n ■ p{n ■ p + uj). The factorization formula can be 
proven in an analogous way to the leading order factorization proof as presented, for example, 
in [38] . The only difference is that we have subleading jet functions instead of a leading jet 
function. 

The second issue, which does not arise at leading order, is the "correct" definition of the 
subleading jet functions and the role of parity and time reversal (PT) symmetry. In general 
the subleading jet function is the discontinuity of TOP of two different combinations of hard- 
collinear fields Oa and O^. Here O can be a hard-collinear quark or a product of hard-collinear 
quark and a hard-collinear gluon, or even a more complicated object. As a result we typically 
have both 

J d^xe"P''{n\T{Ol{0),Obix)}\n) and J d^xe-'P"" {n\T ^^Ol{0),Oaix)^ \n). (50) 
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The subleading jet function(s) should be the discontinuity of the sum of the two terms. In 
practice we find that after we decompose each TOP according to the different color and Lorentz 
structures, we can use translation invariance and the PT symmetry of the strong interaction to 
relate the two TOPs. We can therefore define the subleading jet function as the discontinuity 
of the coefficient of a specific structure in either TOP. 

We will illustrate both issues in the more detailed calculation of the first subleading jet 
function. These details will be suppressed in the derivation of subsequent jet functions. As 
before, the one loop calculation are performed in c? = 4 — 2e dimensions and we implicitly take 
II ^ /xeT^/2(47r)-i/2. 

The following list of the subleading jet functions is not necessarily the complete list of 
possible subleading jet functions. For example, if the QCD currents are matched onto SCET 
beyond tree level we would expect more complicated functions which can depend on more 
than one variable. We wish to emphasize, though, that the list is sufficient to describe all the 
terms of the form of equation (49) at zeroth order in as^i^h) and to all orders in as{iii). In 
particular it includes all the terms of the form hP ■ (8) s° at 0{as). 



5.1 The list of subleading jet functions 

5.1.1 jr. 

This jet function arises from the TOP of the leading order current J*^°^ = e~^"^'''"'^ XT (^S^h^^ 

and the second order current J*^^^ = —e~^"^'''"''^Xr-^n ■ Ahc{S^h}^ . These currents are 
matched at tree level and as a result the Wilson coefficients equal 1. Consequently, the 
hard function, which is simply the product of the Wilson coefficients, equals to 1 also. The 
contribution of this combination of currents to the hadronic tensor is given by 

^Im {B\i J d^xe^(™^)-^T|- (/i5)n-yi,,e^r,X(0), XPj 



^ 2nME 



(51) 



The leading order Lagrangian does not contain any interactions between hard-coUinear and 
soft fields. Since the S-meson states contain only soft particles, we should take the vacuum 
matrix element of the hard-coUinear fields. Prom equation (51) we seem to have two such 
matrix elements, but these are related by the PT symmetry of the strong interaction, as 
explained in appendix A. Consequently, we have 

\hiJn{p')=- j d^xe''P-%n\T{Xk{^\ [Xn-Ahc]^{x)]\Q) 

^ j d^xe-'P-%n\T {[n-AhcX\{Q),Xk{x)]\n), (52) 
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where k, I are color indices. Inserting this definition into (51), we have 



Wij = -Im / d^xe'^'^-""'''^'' 



e 



(27r)^ 



TT 



1 1 



IB\ 



Defining the subleading jet function as 



Up^) = -Im [tJnip'')] , 



TT 



(53) 



(54) 



and using the definition of the leading order shape function [27, 38] 



{B{v)\{hS)^r{S^h)^_\B{v)) 1 



tr r 



1 + 



(55) 



2Mb 2 
we finally find the factorization formula 

Wij^--^ti[r,ir,ifi] [ douUpDsiu), (56) 

which is of the general form of equation (49). 

An explicit one loop calculation of the bare subleading jet function which can also be 
extracted from the sum of equations (43) and (45), gives 



CpOis 



jn{p') = 0{p')^-^ (- + :: + ln^,). 



47r 



1 5 
7 + 4 



fx' 



P^ 



(57) 



5.1.2 jn' 

This jet function arises from the TOP of the leading order current J^°) = e~^"^'''"'^ XT {S^h)_^ 
and the second order current J^^^ = — e"'"*''^'^ X F n ■ Ahc {^^h)^ . We define 



n ■ p 2 



SkiJn'{p')= / d^xe-'P-^iQlTlXkiO), 



X 



1 



/ 



d'^xe-'P-^niT 



n ■ A 



he 



in ■ d 
1 

-in ■ d 



n ■ A 



he 



x) \n) 



X 



{Q),Xk{x)\\Q), (58) 



where A;, I are color indices. We define the subleading jet function jn' as, 

jn'ip^) = -Im [iJuip")] . 

TT 



The contribution of j„/ to the hadronic tensor is 

1 



n ■ p 



tr 



r.^r,i± 

2 ^ 2 



dujjn'ipl) S{u) . 



(59) 



(60) 



An explicit one loop calculation of the bare subleading jet function which can also be 
extracted from the sum of equations (44) and (46), gives 



..(p^) = ^(p^)^-6 0-l + ln| 



(61) 
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5.1.3 and j^^ 



These two jet functions arise when combining two first order currents. Recall from (23) that 
the first order current is 



j(i)^ _e--''--x|^^^, ^ 



in ■ d 



(62) 



Since there are two terms in the first order current, one might naively assume that we will 
need to define a different jet function for each pair of terms. This is not the case, as we explain 
in detail in appendix B. Schematically, the reason is that the inverse derivative acts on all the 
hard-coUinear fields in the first term of (62), which translates into an overall 1/n ■ p factor. 

For all the terms we need to consider the following decomposition of the TOP of hard- 
coUinear fields. 



/^.e-(n|T{ 



4^ T 



(0) 



n ■ p 



91 



2d-2 



^2 ' (i-2 



(63) 



It is clear that g'^ and ie'^^ are the only possible tensors the TOP can depend on. The Dirac 
structure that accompanies each tensor is determined by PT invariance (see appendix A) . We 
now define as usual. 



jf,(p2) = ^Im [i Jlip^)] and j,l(/) = ^Im [i J^,{p')\ . 



(64) 



By an explicit calculation one can show that while jf^ starts ©(ct^), jj^ is non-zero only at 
0{a^g). The one loop result for jf^ is 



(65) 



The contribution of jf^ and j^i to the hadronic tensor is 



n ■ p 
16ml 



tr[r,|ir,^] / dujj^,{pi)S{uj) 



n ■ p 
16ml 



tr [Ti r^- ■ 



du;j^,ipl)S{uj) 



An • p 



tr 



1 + f 



16mh 



tr [r, 7^ r, 7^] - tr [V, 7^5 7^5] 



X 



jiiipD + jtiipl 



S{u). 



(66) 
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5.1.4 jc and Jk 

This jet function arises from the TOP of the leading order current J(°) = e-^"**-^'^ XT {S^h)^_ 

and the second order current J^^) = - e-^^f^'^' XT ^ (¥±^c^±^c) /^t/^) . We find it more 
useful to use the original form of the current as it appeared in [14, 50]. Using the identity 



(67) 



we find that the current can be written as, 

1 {i'P±hc4'±h. 



XT 



in • d rrib ^ in-D nib 

1 [W^ {iD^HcfW + ^W^a,.G^^W] ^^^^^ 



in ■ d 



(68) 



Since {iD^hcY is even under PT symmetry, while G^^ is odd, we need to separate the two 
parts of this current. We define. 



2 1 

w ^X 

—in • d 



(0), Xk{x) \ (69) 



d-2 2' 



75 4;Jg(p')= / dSe-*^-^(l^|T<^Xfc(0), 



m • a 



d^xe-'P-^n\T 



-in • d 



X 



x) I P) 
(0), Xfc(a;)l|f]),(70) 



and the corresponding subleading jet functions. 



Jk(p') = ^Im [i Jk{p^)] and jg(p') = ^Im [i Jg(p')] ■ 
Their contribution to the hadronic tensor is 



(71) 



Wij = tr 

nib 

1 



' 2 ^ 2 



tr[ri^75r,- (^ 



16m5 

An explicit one loop calculation gives. 



dujjK(pl)S(u;) 
?1)75] j dujoipDSiu). 



(72) 
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-2 - + - + In ^ 



(73) 



Notice the interesting fact that at one loop 2jK + in = 0. Similarly 2Jk + Jn = ^, which 
is true even without expanding in e. It is unclear whether this is a one loop "accident" or a 
more general result that holds to all orders in perturbation theory. 

5.1.5 jA and js 

This jet function arises from the TOP of the leading order current J'^^^ = e"*"*''"'^ XT (6*^/1) _^ 
and the second order current J^^) = - g-^"^"^-^ X ^ ^ T id^hc (S^h) ■ We define 

^ ^ ^^^^ 



5w Ja(p'). 



PT symmetry imphes that 



Notice the minus sign in front of the second term. 

We now define the corresponding subleading jet functions, 

jsip") = -Im [i Js{p^)] and ja(p') = -Im [^ Ja(p')] 

TT TT 

Their contribution to the hadronic tensor is 

1 



W, 



16m;, 



tr [r, r, 7^] - tr [r, 7^^75 r, 7^75] 



X J du [js{pl) + {3-d)u{pl)] S{, 
An explicit one loop calculation gives, 



jAip') - eip'f'''' 



Air 



(74) 



(75) 



(76) 



(77) 



(78) 



21 



5.2 Renormalization 



Having defined and calculated all the subleading jet function to one loop order, we are ready 
to discuss their renormalization. We note first that only 3 out of the 8 functions we have 
defined require renormalization at one loop order. These are 



5 
4 



1 + ln 



p2 



(-2) 



1 5 u2 
- + - + In % 
e 4 p2 



In order to renormalize these function we have to introduce a new function 

Jo(/) = %') + 0KO, 



(79) 



^0) 



where only the zeroth order in part of jo(p^) is needed in order to renormalize the subleading 
jet functions at first order in a^. 

It is very tempting to relate this function to the integral over the leading order jet function, 
namely to identify jo(p^) with 



(81) 



(where the lower limit of the integral can be any negative number), since both are equal at 
0[a^g). A very similar function, j{\n.^,^) was defined in [25] 



j(ln4,/^)= r dp'^Jip'^i,). 
A* ^0 



(82) 



In that paper the authors derived the two loop expression for j and its anomalous dimension. 
Since we only use the 0{a^) expression for Jq, we will refrain from identifying it with equation 
(81) or with j of [25]. It seems plausible that all these expression are the same, but in order 
to determine whether they coincide beyond C(a°), a two loop calculation of the subleading 
jet function is needed. Such a calculation is beyond the scope of this paper. 

The most convenient scheme in which to renormalize the subleading jet functions (as well 
as the subleading shape functions) is the DR subtraction scheme [51, 52]. In this scheme the 
Dirac algebra is performed in = 4 dimensions, while loop integrals are evaluated in = 4 — 2e 
dimensions. Choosing this scheme ensures that the renormalized subleading functions are the 
same for B — > I P and the Q7'y — Qjj contribution to B ^ Xgj. 

We renormalize the subleading jet function in the following way. Define the matrix 
^(p^p'^/^) via 



fUp')\ 

in'(P^) 

KMp'')] 



dp'''Z{p\p'-\i,) 



jI^'^p'^) 

J^'%p'^] 



(83) 
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At one loop order we find that in the DR scheme 



I2\ 



/ -4a 1 \ 

-6a 1 
\ 2a 1 y 



(84) 



where a = Cpag/ Ane. Notice that Z {p'^ , p' , ^) is not a square matrix, since we do not 
renormalize jo- If indeed jo is related to the integral over the leading order jet function, then 
Z will have to include other distributions apart from a delta function, for details, see [25]. 

Prom Z(p^,p'^,/x) we can find the renormalization group equations for the subleading jet 
functions. We have at one loop 



dln/i 



in'(P^) 



dp"5{p'-p") 



,2-. CpOis 



At: 



I 8 0\ 

12 
\ -4 y 



( Jo{p")\ 

Up") 
Up") 
\jk{p")/ 



(85) 



The expression for the renormalized subleading jet functions are 

r,2 



Up', /^) = 0{p')^ • 4 ( ^ + In ^ ) + 0{al^ 



P 



Up\^) = o{p')^ 
jK{p',n)-e{p')^ 



(86) 



It is unclear whether with the inclusion of jo the list of the subleading jet functions closes 
under renormalization. In particular one might expect that we need to define a more general 
subleading jet function that depends on more than one variable. For example the discontinuity 
of the Fourier transform of 



{Q\T{Xk{0),n-A'^^{y),Ux)}\n) 



(87) 



where k, I, m, n are color indices. An interesting fact about this type of subleading jet function 
is that at order it only becomes singular in the y ^ Q and y ^ x limits, where it can be 
related to j„. 

Clearly this topic deserves further study, but for our purposes the renormalization via the 
introduction of jo is sufficient at the order in which we are working, namely, as{fii)/mi). 

For phenomenological applications it is convenient to set the scale of the subleading jet 
functions to be the same as the scale of the leading order jet function. This is also the 
convenient scale in which to extract the leading order shape function [29]. At any case the 
resummation of the subleading logs is expected to be a small effect. As a result, the renor- 
malization group eqiiations of the subleading jet function are not expected to be important 
for phenomenological applications. 
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5.3 Results 

We are ready to summarize our results. The subleading jet functions' contribution to the 
hadronic tensor can be written as 



SJF 



du! 



jnipl, jn' ipl, )U) + 3k{pI, IJ')— + joipl, /^) — 

nib ^ ■ P l^b ITT'b 



mt n ■ p nib iT^h 



nil 



T. n n 



n ■ p nib fTT-b 



+ js{pl, fJ'){—- — ]+ JaIpI, iA \ --^^ 



nib rnb 



nib rnb 



S{uj), (88) 



where the traces Ti...Te are^ 



Ti = - tr 
2 



T2 = — tr 
16 



= - tr 



1 



1 + 



f4 = — tr [r,yi75r,-(^-^)75], 



1 



^5 = - tr [r, r, 7^] , = - tr [r, 7^75 T,- 7^75] 



16 

TV = — tr 
^ 16 



Ti yi75 r^- ^75 



(89) 



^Under the assumption that Fj and Tj contain the same number of Dirac's gamma matrices, we can relate 
fi...fe to T1...T4 in equation (12). The relations are Ti = f'i/2, T2 = 71/2 - 2T2, T3 = 2313, = Tk- 
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The subleading jet functions ju, jfi, jn, jn', Jk, jo, js, Ja are defined in section 5.1. The renor- 
mahzed one loop expressions for them are 



3A{P^,^J) 



9{p 



Att 



+ 0{ai) 

e{p') 



2.CFas{ll) 



Air 



0{p') 



47r 



( - 1) + 0{al) 



5 + 41n^ ) +0{a^) 
-6 + 61n^ ) +0{a 



0{p') 



47r 



e{p') 

9{p 



An 



( - 1) + 0{al) 



2 CpOisilj) 



An 



■-]+0{a 



dip 



An 



+ 0{al) 



(90) 



These expressions are in the DR subtraction scheme [51, 52]. As explained in section 5.2 this 
is the most appropriate renormalization scheme for the subleading jet and shape functions. 

The scale dependence in equation (88) cancels against the scale dependence of the sub- 
leading shape functions' contribution [27] 



/SSF 



j duj S{n ■ p + uj) 
v{u;,fi) 



Ta — na^ 



n ■ p 



i2 H -t 1 ' 

rrib 



rrib 



— -ts H = -t 1 



n ■ p 



n ■ p 



-fi H z -'4 



(91) 



n ■ p n ■ p 

For the definition of the subleading shape functions and the traces T1...T4, see section 2. 

We now specialize to the cases of semileptonic and radiative B decays, using the expressions 
for Fj and Vj given in section 2. For the Q77 — Q77 contribution to S — > 7 we need the 
"trace" of the hadronic tensor: W. Its relation to the photon spectrum of S — > 7 is given 
in equation (7). The subleading jet functions' contribution to W is 

4 



W^'^ = Jduj (^Ajk{pI, fi) + AUpI /i) + 2jg{pI fi)^ 

2n • p 



n ■ p 



■Jn'iPoj,lJ^) 



mt 



(92) 



Note that js and Ja do not contribute to W. At the lowest order in as this expression is 



W 



SJF 



— [ dcu — 
J rUb 



1 CpOisip) 



An 



0{pI) 

25 



32 In^ 

Pi 



18 



S{u:) + 0{al), 



(93) 



where in order to simplify the expression, we have used the fact that for this decay mode 
n • p — nib. For completeness we list also the subleading shape functions' contribution 



l^ssF _ / j^^r^ . p _|_ _ ^ S(uj, jji) + siio, jj) — tiio, jji) + uiio, jji) — v(uj, ji) 

rrib 



7iasfu{(^,IJ')-7lOisfvi(^,fi) +0{as) 



(94) 



FoT B ^ X^lP we need the three "form factors" : Wi, W2 and VFcomb = fW^s + W4 + ^#5. 

Wi are defined in equation (8) and their relation to the triple differential decay rate is given 
by equation (9). The subleading jet functions' contribution is 



n ■ p 
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(95) 



At the order in which we are working in, it can approximated by n ■ p/rub- As for W , js and 
jA do not contribute to Wi, W2 or Wcomb- At the lowest order in the last equation is 
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(96) 



If we are using the so called "BLNP" approach [29], i.e. using the definition of y as in equation 
(10), then some of the subleading terms are absorbed into the leading order formula (17). The 
subleading jet functions' contribution to Wi in this case is given by 
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where there is no change to W2 and Wcomb- For completeness we list also the subleading shape 
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functions' contribution 
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(98) 
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6 Summary and Conclusions 

Decay rates of inclusive B decays, namely B X„ / z/ and the Qj^ — Qj^ contribution to 
B ^ Xs'~f, are known to factorize in the end point region at the leading order in Aqcd/""^?) 
into a product of a hard function and a universal leading order jet function convoluted with a 
universal leading order shape function. The hard and jet function are calculable in perturba- 
tion theory while the shape function is a non perturbative object. Recently the hard function 
for semileptonic decays was calculated at O(a^) [15, 16, 17, 18]. Together with the already 
known two loop calculation of the jet function [25], this will allow us to reach a full 0{a'^) 
accuracy, at leading power, in describing the decay rates in general, and in extracting \ Vub\ in 
particular. 

Beyond leading order in Aqcd/"^6, one would expect the decay rate to factorize into sums of 
products of subleading hard functions and subleading jet functions convoluted with subleading 
shape functions. Of these power suppressed corrections, only the subleading shape functions 
were known. In this paper we have analyzed the subleading jet functions' contribution. These 
arise first at order 0{as) and appear in the partial rate convoluted with the leading order 
shape function. 

First, we have argued that at order AqcD/^Tib only subleading jet and shape functions 
contribute. Subleading hard function can only appear when multiplied with subleading jet or 
shape functions. We have demonstrated this explicitly at one loop by showing that the 0{as) 
corrections which are Aqco/^b suppressed in the end point region, arise from two momentum 
regions: a hard-coUinear region and a soft region. We have then shown that the soft region is 
accounted for in the parton model by the one loop matrix elements of the known "tree level" 
subleading shape functions. The hard-collinear region is accounted for in the parton model 
by the time ordered products of subleading SCET currents that apart from the heavy quark 
itself, do not include any soft fields or soft covariant derivatives. 

In the main section of this paper, section 5, we have defined to all orders in cts for the 
case of a tree level hard function, the 8 subleading jet functions that can contribute to partial 
rates of inclusive B decays, and calculated their one loop expressions. After a short discussion 
of the renormalization of the subleading jet functions, we specialized to the phenomenologically 
interesting cases of S — > Xu I P and the Q77 — Q77 contribution to B ^ Xsj and presented 
explicit expressions for these decay modes. The main results of the paper are collected in 
section 5.3. 

We can now summarize the factorization formula for inclusive B decays in the following 
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way 

known new 




where the label "new" refers to the new results of this paper. The label "known" refers to terms 
in the factorization formula for which we have explicit expressions for all of their perturbative 
components. Thus H is the leading order hard function, J is the leading order jet function, 
both known at 0{af), Jq is the C(q;°) part of the leading order jet function, h — 1 + 0{as), 
and jk are given in section 5.3. The rest of the l/rub suppressed terms for which we do not 
have such explicit expressions can be found in [26] and [28]. 

While the subleading jet functions' contribution is both as and l/rrib suppressed in the end 
point region, their contribution becomes more important as one moves out of the this region, 
since the l/m^ suppression is reduced as one is integrating over larger and larger portions of 
phase space. The one-loop subleading shape functions contribution, although formally as/mj 
suppressed in the end point region, is expected to become even more power suppressed outside 
of the end point region. Furthermore, the kinematical area outside of the end point region is 
becoming more important with the constant improvement of experiments and the relaxation 
of kinematical cuts. Together with the recent calculation of the leading order hard function for 
semileptonic decays at 0{al), this work takes another step towards a more precise description 
of inclusive B decays in the end point region. 

Although we have not presented any kind of numerical analysis, the implementation of the 
subleading jet function analysis with the framework of "BLNP" [29] is relatively easy One 
needs to replace equation (15) by (94) and (16) by equations (96) and (97). At the same 
time one needs to modify the treatment of the subleading shape function in [29]. Only the 
combined numerical analysis would be meaningful. Such a study is left for future work. 

Another issue that deserves further consideration is the renormalization of the subleading 
jet functions. In particular, it would be desirable to have a complete basis of subleading jet 
functions. As in the case of the subleading shape functions [49] , it seems that further study is 
needed if we wish to understand the renormalization and mixing of these non local operators. 

Finally, subleading jet functions arise also outside of flavor physics, for example, in the 
X ^ 1 region of deep inelastic scattering. More spcciflcally, the subleading jet function jf^ 
appear in the factorization formula for the longitudinal structure function [33, 34, 35, 36]^. A 
more detailed study of the subleading jet functions' contributions to the x ^ 1 region of deep 
inelastic scattering is left for future work. 
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by the United States-Israel Bi-national Science Foundation grant 2006280. 

'^Note though that the subleading jet function in [33, 34] and [35] differ from our definition and from each 
other. 
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A Consequences of PT symmetry 



In this section we will prove our claim that PT symmetry and translation invariance allows 
us to relate the TOP of two different hard-coUinear operators. Define 

Ta,^ j d^xe-'^^ {n\T {Ol{Q),0,{x)]\n) ^ 

d^xe-'^'^[^(x°)(O|O6(x)Ol(0)|O)±^(-a;°)(O|Oi(0)O6(x)|O) 



d^xe-'P'' 



^(x°)(O|O„(x)Oj(0)|O) ± e{-x'){n\Ol{0)Oa{x)\n)] (99) 



We would like to relate T^b and T^a- Using translation invariance, and the PT invariance of 
the strong interactions and of the vacuum, we have 
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In order to relate T„h and T/,„ we need to relate 
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As our first example consider In section 5 we defined. 



5) hlJnip') 
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(102) 

where A;, / are color indices and a, h are spinor indices. Consider now the other TOP: 
Tother^ j d^xe-'P-%n\T{[n-Ahc'Aai{Q),%k{x)]\Q)^ 



I 



d^xe-'P-"" 



9{x''){Q\Xbk{x) [n-AhcXli (0)|fi) - e(-x°)(Q| [n ■ AhcXli {0)%k{x)\n) 

(103) 
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Using translation invariance and the PT symmetry we have, in the Weyl representation of 
Dirac 7 matrices, 



other 



9{x'){n\ [X^'^''n-Anc],,{x) [7^7' X],, (0)|1]) 
-^(-a;°)(l]|[7'7'X],,(0)[X7S'n-^,e]„,(a;)|l^) 1. (104) 



Comparing (102) and (104) and using y 1 1^ T 1 — il^) we find that 



ha 



(105) 



ab 



which proves equation (52). 

As a second example, consider jf^ and J-^. Here we use the same transformations to prove 
the decomposition the TOP into two jet functions. In section 5.1.3 we defined 
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Using translation invariance and the PT symmetry we have 
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(107) 



which proves that is accompanied by (?^/2), while ie^^ is accompanied by (^75/2). In a 
similar manner we can use PT symmetry to analyze the rest of the TOPs. 



B T{ J(i), J(i)} subleading jet function(s) 

In section 5.1.3 we argued that the TOP of the first order SCET current with itself, gives 
rise to only one subleading jet function. The reason is that the inverse derivative is acting 
on all the hard-coUinear fields. On the other hand, we were forced to defined two different 
jet function j,„ and j„' since the inverse derivative was only acting on the hard-collinear gluon 
and not on the hard-collinear quark. In this appendix we give a rigorous proof of these facts. 
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We consider the TOP of two currents: Ji{x) and j^J2{x). Define 
T2i = ij d^xe-'P-^'iniT yi{0), Ji{x)^\fl) 

T2>i = i [ d'^xe-'P''{n\T \jI{0) Ji{x)\ |Q). (108) 

J I —in ■ d ) 



We would like to prove that 



^ Im (t^2' + T2>A = -3^- Im ( + ) . (109) 



(110) 



TT \ / n ■ pn 

Consider T12' first. Inserting a complete set of states |r) we can write it as 

T,2'^i [ d'xe-'P-^ \ e{x')J2{^\^M^)\r){r\Jl{0m) 

^ r 

±e{-x') ^ (O|4(0)|r)(r| ^ J2{xm 

r 

Using translation invariance we find 

Ti2'=i I d^xe-'P-"" ^(a;°)y 3^(1^1 J2(0)|r)(r|4(0)|l])e-^^^ 

E (Z^ (O|4(0)|r)(r| J2(0)|O)e-^j (111) 
(we use r to denote both the state and its momentum). Using the identity 

e(x^)^^ du- — , 

and integrating over x and lu, we obtain 
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Repeating the same procedure for T2'i we find 
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Since that by definition {a\Ji\b)* = {b\jj\a), we find that 
Tw + T2'i = ^ 2 Re 

r 

±^2 5^2Re 



n ■ r 



(Q|J2(0)|r)(r|4(0)|f})^ 
{n\j\mr){r\Um) 



(27r)353(p + f) 
-Po - ro + ie 



(— n • r) 



(27r)3(5^(p - 
Po - ro + ie 



(114) 



Using the identity Im [1/(m + ie)] = —t:6{u), we have 



-Im(Ti2'+T2a J = V2Re[(l]|J2(0)|r)(r|4(0)|fi)l (27r)='5^(p + r) 
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(115) 



We should note that throughout the above derivation we have assumed that the n component 
of the total momentum of the state is not zero and therefore its inverse is defined. This 
assumption follows from the fact that by definition the hard-coUinear states have a "large" n 
component of momentum. In the same way we can prove that 



-Im(i / d'^xe-'P-^(n\T \4(0) , J2(x)\ IQ) 

n \ J { _ifi. d in-d J 
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(116) 



Using these identities we obtain the results of section 5.1.3. 
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